Acoustic Analogue of Bloch Oscillations and Resonant Zener Tunneling in Ultrasonic 

Superlattices 
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We demonstrate the existence of Bloch oscillations of acoustic fields in sound propagation through 
a superlattice of water cavities and layers of methyl methacrylate. To obtain the acoustic equivalent 
of a Wannier-Stark ladder, we have employed a set of cavities with different thicknesses. Bloch 
oscillations are observed as time-resolved oscillations of transmission in a direct analogy to electronic 
Bloch oscillations in biased semiconductor superlattices. Moreover, for a particular gradient of 
cavity thicknesses, an overlap of two acoustic minibands occurs, which results in resonant Zener-like 
transmission enhancement. 

PACS numbers: 43.20.+g, 43.40.+S, 46.40.Cd 
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Bloch oscillations (BO) and Zener tunneling arc fun- 
damental transport effects appearing when electrons in 
a periodic potential are accelerated by an external DC 
electric field [l], Both phenomena have been impres- 
sively demonstrated in a number of experiments after 
the advent of semiconductor superlattices (SSL)[3(. The 
frequency domain counterpart of BO is the equidistant 
Wannier-Stark ladders (WSL) of the electronic states in 
a biased superlattice, leading to resonances of the den- 
sity of states which were observed in optical spectra^. 
In time-resolved optical experiments, BO were first ob- 
served as oscillations of electron wave packets in biased 
SSL [H, [H, 0, @, and later as a periodic motion of 
ensembles of ultracold atoms [lot 11 1 and Bose-Einstein 
condensates [13, Ha ]. The related high-field phenomena, 
the nonresonant Zener tunneling between neighboring 
minibands and resonant Zener- like tunneling between the 
anticrossing Wannier-Stark states of neighboring mini- 
bands, were also observed in SSL 14. 1511. 



In general, the Bragg reflection can cause BO of a wave 
of any nature (electronic, optical, acoustical or matter 
wave) in a lattice with a weak linear gradient of the lat- 
tice potential, which can be caused in turn by an external 
field or perturbation of any nature (electric, magnetic, 
acceleration or gravitation field), see, e.g., [lq |. Opti- 
cal BO and Zener tunneling of light waves have been 
recently observed in time-resolved experiments in 2D(l7j 
and lD[l8| optical superlattices with refractive index gra- 
dient along the growth direction. Phonon BO and Raman 
spectra of phononic WSL states were described in semi- 
conductor multilayer solid structure based on acoustic- 
phonon cavities with different thicknesses [19]. More re- 
cently, WSL have been also observed in one-dimensional 
elastic systems (2fj|. 

In this Letter we predict analytically and confirm ex- 
perimentally the existence of acoustic Bloch oscillations 
in two-component ultrasonic superlattice made of layers 
of Methyl Methacrylate (Plexiglas) and water cavities. 
The system of Plexiglas' layers, all having the same cho- 



sen thickness, separated by water cavities with a given 
gradient of thicknesses, is equivalent to a set of acoustic- 
wave cavities with a gradient of local resonant frequen- 
cies. In this system the water cavities act as resonant 
cavities, the resonant acoustic modes being well confined, 
which results in well-pronounced acoustic WSL in the 
transmission spectra through the structure. This in turn 
will result in time-resolved Bloch oscillations in the trans- 
mission (and reflection) spectrum for an incident acous- 
tic pulse with the proper spectral position and width. 
We will also show that for a particular gradient of water 
layer thicknesses in the superlattice, an overlap of two 
acoustic bands occurs which results in resonant enhance- 
ment of the phonon transmission through the superlattice 
due to resonant Zener-like effect of spatially-overlapping 
phononic states belonging to neighboring minibands. 

We start with the evolution equation for the Bloch 
wavevector kz (along the axis of the periodic system) : 



• dkz duo ■ du 
Z = IfadZ dZ* 



(1) 



where u is phonon frequency and dcu/dkz = vz = Z 
is the group velocity of the Bloch phonon wave packet. 
Equation (1) presents a generalization for phonons of 
the "acceleration theorem" for the evolution of elec- 
tron Bloch wavevector k in an external electric field E, 
fik = eE = F, with duj/dZ playing the role of the exter- 
nal "wave" force Fz/h. 

The dispersion equation for the longitudinal Bloch 
phonons (sound waves) propagating along the z— axis 
of two-component superlattice is the following, see, e.g., 
Ref. [5l|: 



cos kzd — cos( ) cos( ) 
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where d,A,w, PA,w and ca,w are the thicknesses, densi- 
ties, and sound velocities in the layers of elastic materials 
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FIG. 1: (Color online) Transfer-matrix calculation of the in- 
tensity distribution (log|P(,z)| 2 ) through a sonic superlattice 
made of eight coupled water cavities, (a) Flat band case, 
A(l/dw) =0%, where two minibands (MB1 and MB2) are 
clearly seen, (b) Tilted band case, /S.(l/dw) =6%, where 
Wannier-Stark like ladders (WSL1 and WSL2) are also shown. 
The structure of coupled water cavities is schematically shown 
above each panel; dark regions represent Plexiglass layers. 
The right panels depict the transmission spectra across the 
total structure. 



A and W, d = (La + dw is the period of the superlattice. 
No mode conversion (from longitudinal to transversal) 
has been considered in this problem since the propagating 
sound wave is assumed that impinges the elastic layers 
normally to their surfaces. 

In order to separate the Fabry-Perot resonances of the 
transmission through the weakly-coupled elastic A and 
W layers, let us also assume that there is big acoustic 
mismatch between layers A and W; i.e., pwcw "C Paca- 
The corresponding unperturbed systems have their reso- 
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n and m integers. ^From Eq. (2) it is easy to get a nar- 
row isolated tight-binding phonon miniband centered at 
the first (m = 1) Fabry-Perot resonance in W layers by 
imposing the condition cw/dw < CA/d a - For example, 
for cw/dw = ca/^a one has (in the reduced Brillouin 
zone —7:<kzd<Tr): 



c w 
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dw dw PACA 



(3) 



Now, according to Eqs. (1) and (3), for a constant 
phonon driving force Fz/h— duj/dZ the phonon group 
velocity dto/dkz will oscillate with the Bloch frequency 
u}B=Fzd/h=ddcu/dZ, in a direct analogy to electronic 
Bloch oscillations in a biased SSL. 

To demonstrate the Wannier-Stark like ladders and the 



corresponding Bloch oscillations predicted by the ana- 
lytical model above, we have considered a superlattice 
made by alternating fluid layers and layers of elastic ma- 
terial. In particular, experimental measurements have 
been performed on a set of multilayers (ML) samples con- 
sisting of 8 coupled water cavities enclosed by 9 lay- 
ers of Plexiglas. Transfer-matrix (TM) calculations [22| 
of the sound intensity transmitted through this multi- 
layer as a function of frequency have been performed by 
using the following data. As parameters we have taken 
Pa =1.19 g/cm 3 , ca =2.65x 10 5 cm/s and pw =1 g/cm 3 , 
cw =1.48xl0 5 cm/s for the density and sound velocity in 
Plexiglas (layers A) and water layers, respectively. Vis- 
cosity effects are physically irrelevant to the studied phe- 
nomena and are not reported here. 

First, a perfect superlattice consisting of alternating 
layers with equal thickness, dA = dw =0.16 cm, has been 
studied by TM simulations. The results for this case (see 
upper panels in Fig. 1) show that, though the contrast of 
acoustic impedances between Plexiglas and water is low, 
a group of modes MB1 is fairly well confined in the water 
cavities. The localization of these modes is visible as a 
set of elongated hot spots surrounded by dark regions 
representing the bandgaps. The group is centered at the 
first (m =1) Fabry-Perot resonance of water cavities and 
corresponds to a discrete analog of a miniband [see Eq. 
(3)]in our eight-cavity structure. The liner midfrequency 
of MB1 obtained from TM calculation is uf M =446 kHz, 
which agrees with the first Fabry-Perot cavity resonance 
in the water cavities; i.e., v\ ~ (l/2)(cw/dw) =463 kHz. 
The localization of the states belonging to the group MB2 
is not as clear as in MB1 because MB2 contains states 
that are also localized in the Plexiglass layers. 

In our acoustic ML structures the frequency of states 
localized in the water cavities depends approximately lin- 
early on the driving force, duj/dZ, with the slope given 
by mdwduj/dZ (m =1, 2, 3, ...), where dcu/dZ = 
■KCw-§z{^/dw)- For example, for the case m =1, a lin- 
ear variation of frequency is obtained by introducing a 
constant variation of the inverse of the thicknesses of the 
water cavities 1/dw [see Eq. (3)]: 

6ui = Ld(ze) - w(ze-i) = nc w S(-y—) (4) 

dw 

where 6(1/ dw) = (l/dw e — 1/dwe-i) = const, the 
subindex I defines the ordering of cavities in the struc- 
ture, I = 1,2,... 8. From here onwards, the results 
will be given in terms of the dimensionless parame- 
ter A(l/d w ) = [(l/d We ) - (l/dw<_i)]/(l/dwi), where 
dwi =0.16 cm is the thickness of the first, £ =1, cavity. 
Figure 1(b) shows the pressure intensity distribution for 
A(l/dw) =6%. Notice that the group MB1 is linearly 
tilted when a linear gradient is imposed over the cav- 
ity thicknesses. It becomes in a discrete sequence of fre- 
quency levels, WSL1, strongly localized within one of the 
water cavities and one therefore obtains a well-defined 
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FIG. 2: (Color online) (a) Amplitude of transmission calcu- 
lated by transfer matrix (TM) for the multilayer described in 
Fig. 1(b). The acoustical analogue of the electronic Wannier- 
Stark ladder is seen as a series of almost equidistant transmis- 
sion peaks, (b) The corresponding experimental spectrum ob- 
tained from time-domain measurements in transmission, (c) 
Spacing between peaks (&vb) in the corresponding WSL as a 
function of the gradient A(l/dw)- The gray dotted lines are 
guides for the eye and allow to distinguish the linear WSL 
regime. 



spacing 8u>bi between levels. This set of equally spaced 
levels are the acoustical equivalent of the Wannier-Stark 
levels. The spacing between levels being approximately 
determined by Slubi — u(zi)A(l/dw)- The discussion 
about the levels in MB2 is more complex, but a set WSL2 
with level spacing Slob2 can also be defined. In what fol- 
lows we will focus our discussion in the set WSL1. 

Therefore, for the gradient considered (6%) the fre- 
quency separation (between states in WSL1) predicted 
by the simple analytical model previuosly described is 
Svbi «27.8 kHz, which is close to that found between 
peaks in Fig. Ufa), which shows the transmission spec- 
trum calculated by TM for the case of a gaussian beam 
f(oj) centered at 500 kHz crossing the gradient superlat- 
tice. A system having such discrete sequence of frequency 
levels with level spacing u>b is the acoustical equivalent 
of the electronic WSL, and is expected to exhibit acous- 
tical Bloch oscillations of period Tg — 2tt/ujb as it has 
been experimentally demonstrated. 

Transmission measurements have been performed in 
a water tank by a simple experimental set up con- 
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FIG. 3: Temporal response of the system for various values 
of the gradient. The top panel shows the unperturbed probe 
pulse without superlattice sample. The observed oscillations 
are the acoustical counterpart of the electronic Bloch oscil- 
lations. The period, Tb, measured are 41.5 fis, 32.7 /js, and 
29.3 /is for gradients 6%, 8%, and 9%, respectively. 



sisting of two transducers, an ultrasonic square wave 
pulser/receiver and a digital oscilloscope connected to a 
computer to store the data. In brief, the emitter trans- 
ducer, which is put near (lees than 1 cm) the sample's 
surface, is employed to excite a longitudinal wave that 
crosses the structure and is recorded by a the receiver 
transducer, which is also placed near the last layer of 
Plexiglas. The central frequency of both transducers is 
500 kHz, their bandwidth being approximately 50%. The 
temporal resolution is 0.4 /is. Effects of mode conversion 
in this system can be considered as negligible since the 
excitation and detection involve only vibrations normal 
to the interfaces and the elastic layers are very thin. 

Figure 2(b) shows, as a typical result, the experimen- 
tal spectrum for the amplitude of states WSL1 corre- 
sponding to a gradient of 6%. The blue arrows indicate 
the position of the peaks determined by TM simulation 
(see Fig. 2a). The frequency spacing between peaks is 
Su B x f «24.3 kHz, which agrees with the value calculated 
by TM 5u^f «24.8 kHz, and with the one obtained by 
the analytical model. Also, an overall good agreement is 
obtained between theoretical and experimental spectra 
in the full range of frequencies of interest. 

The comparison between measured frequency spacings 
6vb as a function of the gradient with the ones calculated 
by TM is shown in Fig. [2fc). It is remarkable the very 
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FIG. 4: (Color online) (a) Transfer matrix (TM) calculation 
of the transmission spectra of a tilted superlattice around 
the value of gradient, A(l /dw)), where the first anticross- 
ing of the acoustical WSL and hence Zener-like resonant ef- 
fect occurs. The amplitude is plotted as a function of the 
reduced frequency, where wo is the central frequency. The in- 
set shows the strength (see text) of the resonant peak arising 
from the interaction between the two states located in WSL1 
and WSL2, respectively [see Fig. 1(b)]. (b) The transmission 
spectra measured by time-resolved transmission experiments. 
The inset shows a comparison between the behavior obtained 
by TM simulations (full lines) and experiments (symbols with 
error bars). 



good agreement between theoretical predictions and mea- 
sured values. It must be pointed out that the true WSL 
regime is characterized by a linear dependence of 5vb as 
a function of gradient; i.e., Svb — * when A.(l/dw) —> 0. 
To check this behavior, two straight dotted lines are plot- 
ted in Fig. [2^c) to allow the distinction between the non- 
linear Fabry-Perot regime. It can be concluded from Fig. 
I2tc) that the linear regime is shorter for WSL2 than for 
WSL1. 

The time resolved transmission experiments are per- 
formed by sending the short pulse depicted in the upper 
panel of Fig. [3J The Bloch oscillations corresponding 
to previously described WSL1 arc also shown in Fig. 3 
for various values of the gradient. The oscillation period 
Tb experimentally observed decreases while increasing 
A(l/dw) as expected. Our previous analysis [see Fig. 
2(c)] concluded that above 6% gradient a true WSL is 
formed in our sample. ^From the time-resolved data we 



can observed that indeed oscillations occur in transmis- 
sion. More than five period oscillations are observed in 
the transmitted intensity with a period that decreases 
as A(l/d\v) increases. In addition, as the gradient in- 
creases the intensity transmitted decreases, which can be 
understood from the increases tilt of the bandgap (see 
also Fig. 1). The experimental data are in very good 
agreement with the calculated dependence, which is not 
shown here but can be directly concluded from Fig. [2fc) 
where 8vb = 1/Tb 

To examine the Zener-like resonant effect, we per- 
formed TM calculations for gradient larger enough to ob- 
serve that one level from WSL1 and another from WSL2 
lead to resonant transmission through the superlattice. 
Figure 4(a) depicts the results showing that an important 
feature of the interaction between levels is its anticross- 
ing behavior. The strength of the interaction, which is 
defined as the integral of the transmission peak, is shown 
as a inset in Fig. [4ja) and its maximum value is achieved 
when the separation between peaks is minimum; i.e. at 
the gradient of 9.93 %. Figure 4(b) shows the transmis- 
sion spectra obtained from time-resolved transmission ex- 
periments confirming the TM predictions. It is possible 
to observe how the two states start to overlap and the 
saddle like curvature transforms into a sharp resonance, 
where the separation between levels is difficult to appre- 
ciate mainly due to the precision of our experimental set 
up. The inset in Fig. 4(b) plots the comparison with the 
TM simulations, the error bars showing that the main 
source of error comes from the precision in the definition 
of gradients, which comes from the mechanical inaccura- 
cies associated to the construction of water cavities with 
correct thicknesses. 

In summary, acoustic Bloch oscillations and the acous- 
tic analogue of the resonant Zener-tunneling have been 
demonstrated in composites consisting of water cavities 
surrounded by layers of Plexiglas (methyl methacrylate) . 
This is the most simple system used to date to observe 
acoustic Bloch oscillations showing that the phenomenon 
is indeed universal and very robust. We hope that this 
demonstration stimulates future work to create acoustic 
devices based on these phenomena. 
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